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ABSTRACT 

The optimum shape of a blunt forebody attached to a symmetric wedge or cone is de- 
termined. The length of the forebody, its semi-thickness or base radius, the nose radius 
and the radius of the fillet joining the forebody to the wedge or cone are specified. The 
optimum shape is composed of simple curves. Thus experimental models can be built read- 
ily to investigate the utilization of aerodynamic heating for boundary layer control. The 
optimum shape based on the modified Newtonian theory can also serve as the preliminary 
shape for the numerical solution of the optimum shape using the governing equations for a 
compressible inviscid or viscous flow. 


1 Research was supported by the National Aeronautics and Space Administration under NASA Contract 
No. NAS1-18605 while the second author was in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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1. INTRODUCTION 


Optimum shapes of two dimensional or axi-symmetric bodies in hypersonic flight were 
studied in the early stage of hypersonic flow research. The results were summarized in the 
book [1] by Hayes and Probstein. The analyses were based on the Newtonian- Busemann 
(N-B) theory which was justified as the leading term of an asymptotic solution at large 
upstream Mach number. Shapes of minimum wave drag with given body length, thickness 
and/or volume were determined. The minimum wave drag is achieved when the product 
Pi cos 9 \ is a maximum where 6 \ is the inclination of the surface at the rear end with respect 
to the free stream while Pi is an integral along the surface from the front to the rear and 
depends on the shape. With the assumption that 61 and Pi are independent of each other, 
the “absolute” optimum shape which maximize the integral Pi has to fulfill also the end 
condition 

= 0. (1) 

This end condition cannot be met without a discontinuity in slope. Engineering designs 
using an absolute optimum forebody followed by a body confined within a free surface were 
proposed to achieve a “proper” shaped body with a wave drag somewhat higher than that 
of the optimum shape. We also note that the optimum forebody has a sharp leading edge 
or a conical nose [1]. 

It is well known that the rate of heat transfer to a “cold” wall at the stagnation point is 
inversely proportional to the square root of the nose radius [2-4]. An optimum shape with a 
sharp or pointed forebody is not suitable for hypersonic flight. To reduce the aerodynamic 
heating to an allowable limit, a body in hypersonic flight has to have a finite radius of 
curvature, i2o, at the stagnation point. The minimum radius depends of course on the flight 
Mach number the altitude and the physical properties of the nose with or without a 
heat shielding device. An optimum forebody has to have a nose with a prescribed radius of 
curvature greater than the minimum. 

In recent studies on boundary layer control by surface heating or cooling [5-8]. It was 
found that a boundary layer which is heated upstream and/or cooled downstream is more 
stable than that without thermal control. In a hypersonic flow, the aerodynamic heating 
near the stagnation point provides a nature source for heating the boundary layer. The 
amount of heat acquired by the boundary layer is proportional to , where j = 0 and 1 
for two dimensional and axi-symmetric flow respectively. From the point of view of boundary 
layer control, it is also desirable to have a nose with a prescribed radius of curvature Rq and 
then find the optimum shape. 

It is well known that the modified Newtonian theory [1], although an empirical theory, 
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has been employed in engineering designs because its prediction of surface pressure is in 
much better agreement with the experimental data than the N-B theory. We shall use the 
modified Newtonian theory to compute the wave drag and find the optimum shape in this 
report. To avoid flow separation or shock boundary layer interaction, we require that the 
optimum contour should be smooth, that is, it should match with the circular nose in the 
front and the wedge or cone in the rear without a discontinuity in slope. As we shall see 
in the next section, the optimum shape under the modify Newtonian theory does not call 
for the rear end condition (1). However, the slope of the optimum forebody in general will 
not match with that of the wedge or cone. To remove the discontinuity in slope we assume 
that the forebody is joined smoothly to the wedge or cone via a fillet of radius R\. Thus 
we complete the specifications of the free front and rear end conditions for the optimum 
contour. 


2. THE OPTIMUM SHAPE 


The modified Newtonian theory is the direct application of the Newtonian theory along 
the surface with a correction factor such that the correct stagnation pressure p a is obtained. 
The surface pressure p is related to the inclination of the surface with the upstream flow, 0, 
by 

p = p a sin 2 9. (2) 

The wave drag of the forebody is 


D = 2p a f sin 2 d (7r y) 3 dy, 
Jo 


(3) 


where h is the semi-thickness or base radius. With dx = dy cot (9, the length of the forebody 
l yields a constraint, 

i = y cot 6 dy. (4) 

Figure 1 shows the front and the rear end of the forebody located at x = — l and x = 0 
respectively. The centers for the circular nose and the fillet are located at Co(—l -f- i2o,0) 
and Ci(— /2sin/3, h + R cos /?), where /3 denotes the half wedge angle or cone angle. Here we 
have assumed that the nose is convex with R 0 > 0 while the fillet can be concave for R x > 0 
or convex for R x > 0. 

The parametric equations for the circular nose and fillet are 

x = (- 1 + Rq ) - Rq sin 9 n , 

= R 0 cos9 n , 


y 


2 


(5) 

( 6 ) 



for 0 < y < yo and 

x = — Ri sin/3 + Ri sm0f (7) 

y = (A. + R x cos/3) — Ri cos Of (8) 

for h > y > y\. The shape function 9 n (y) along the circular nose and Of(y) along the fillet 
are defined by (6) and (8). The shape function of the contour, 9(y) for yo > y > yi, and 
its two free ends, yo and yi, shall be defined by the necessary conditions for minimum wave 
drag subject to the constraint (4). We introduce the Lagrange multiplier A and then seek to 
minimize 

$(ei, £2, e) = -\l= f dy (sin 2 9 n - A cot 9 n ) 

Zp a ir J Jo 

+ f dy (sin 2 9 — A cot 9) + / dy (sin 2 9j — A cot 9f ), (9) 

>'{*> Jvi 

where ex, 62 and 677(3/ ) denote the variations of 3/0, yi, and 9(y), i.e., 

Vo = Vo + Co, Vi = yi + Cl and 


y = Kv) + € v(y) for Vo < y < y\- 

(10) 

The necessary conditions for a minimum $ are 

8{yt ) - 0 n{y o ), 

(11) 

Kv\) = 9 f(vt ) 

(12) 

and yi sin 3 6 — constant hC. 

(13) 


Condition (4) is fulfilled on account of (5) and (6). Equations (11) and (12) are the natural 
end conditions, requiring that the optimum contour joins smoothly with the circular nose at 
yo and the fillet at y\. Equation (13) says that 9 remains constant for yo < y < yi for two 
dimensional cases and sin 8 ~ 3/ _ 1/1 3 for axi-symmetric cases. 

The optimum forebody, defined by (5-8), (11), (12) and (13), depends on four parameters 
Ro, Ri,l and /? with 

Ro = Ro/h, i?x = Rjh and I=l/h. (14) 

In order to study the ranges of these four parameters and their effects on the wave drag, 
we introduce a reference forebody as shown in Fig. 2. It is composed of the extension of the 
wedge (or conical) surface matched with the circular (or spherical) nose of radius Ro. Hence, 
we have the shape function 

9*(y) =13 for h > y > yi = yo = Rocos /3 (15) 

and 9*(y) = cos _1 (y/i2o) for y£ > y > 0. (16) 
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The scaled length of the forebody is 


F = cos/? + Ro (1 — cos /?). 


(17) 


This configuration implies that Rq cos /? < 1 . This is fulfilled in real problems for which 
Ro < 1. Its scaled wave drag is 

= IT 7 [1 " (fl ° coa ' 5),+ ’ 1 + " Trf 1 ' (18) 

Also shown in Fig. 2 is an optimum two dimensional forebody with the same /?, Rq and 
h and a convex fillet, (Ri < 0). The optimum shape is convex with 


0(0) = 7t/2, 6(h) = /? and 8'(y ) < 0, 


(19) 


for 0 < y < h. In particular, we note that at the point of tangency with the circular nose, 

8{y 0 ) > /? hence y 0 < y* (20) 

and that 

6(y) > cos -1 (j//Ro), (21) 

for y 0 < y < y*. It follows from (19), (20) and (21) that 

9{y) > 0*(y) for 0 < y < h, (22) 

and from (3) and (4) that, 

D > D* and l<l*. (23) 

Thus we conclude: 

I, The wave drag of an optimum forebody with a convex fillet , R\ < 0, is greater 
than the wave drag D* of the corresponding reference forebody, i.e., having the 
same /?, Ro and h, while the length of the former is shorter than the length l* of 
the latter. 

The reference forebody can be considered as the limiting optimum forebody as Ri — > — oo 
while y x — > y 0 - 

From hereon, we consider only the optimum forebody with a concave fillet, i.e., 

Rx > 0 (24) 

and study the drag reduction relative to the drag D * of the corresponding reference forebody 
as Ri decreases and l increases. We shall treat the two dimensional and the axi-symmetric 
problems separately in the following two sections. In each section we determine the following: 
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1. The optimum contour and the points of the tangency of the contour with the circular 
nose and the convex fillet for a given set of Ro, R\, I and /?. 

2. The upper and lower bounds of the scaled length T as functions of Ro, Ri and /3. 

3. The dependence of the wave drag on those parameters. 

3. TWO DIMENSIONAL OPTIMUM FOREBODIES 

For j = 0, Eq. (13) for the optimum contour yields 

6 = constant a for y 0 < y < y\. (25) 

The optimum forebody is composed of the circular nose, for 0 < y < yo, and the concave 
circular fillet, for y\ < y < h, joined together by their common tangent ToTi, where To 
and Ti denote the points of tangency at y = yo and y = t/i, respectively, (see Fig. 2). 
Instead of solving for a,y 0 and y x from Eqs. (5-8), (11), (12) and (25), we determine them 
directly by making use of elementary geometry. Let CoC% be the straight line parallel to 
T 0 Ti intersecting the radial line C\T\ at C 2- Then points C0C2 and C\ are the vertices of a 
right triangle with 

\C 1 C 2 \ = R <) + R 1 . (26) 

Let ip and r denote LCiCoC\ and LOCqC\ respectively. We obtain from Fig. 2 the following: 

ip = sin 1 [(I?i + -T2o )] / 1 CqC\ | , 

r = sin -'[(h + Ri coBfi)]/\C 0 Ci\ t (27) 

/S>a = T — V>> 0, 

with ICo^f/h 2 = {l-Ro-R 1 sin Pf + (1 + Ri cos/3) 2 . 

Note that the angles ^>,r and a are functions of 1,Rq,Ri and /? and ip < 0 in Fig. 2 since 
—Ri > Ro > 0. The solutions are real so long as |C7 0 C7i | > Ri + Ro- 

When l is equal to I* of the reference forebody, Eq. (17), we have a = (3 and y 1 = h for 
all Ri > 0, keeping in mind that RoCos/3 < 1. For an optimum contour we have a < (3, 
yi < h and hence l > I*. 

If the ordinate of the center C\ of the fillet is greater than or equal to R\ + R2, i.e., 

1 + Ri cos /3 > Ri + Ro, (28) 

it is evident from (27) that |C'oC'i| > Ri + R 2 and a > 0 and that as l increases from J* to 
00 , a decreases from j3 to 0. 
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On the other hand if (28) is not true, we have an upper bound for l defined by the 
condition \CiCo\ — Rx + Rq. The upper bound is 

lm = R\ sin/3 + Rq — {(i2i + Rq) 2 — (1 + Rx cos/?)} 1 ^ 2 . (29) 

When the length reaches its maximum, we have 6 0 = 6\ = 9 m = cos -1 [(l + Rx cos /3)/(Rx + 

*»)]• 

As the length l increases from l* to l m , decreases from /3 to 9 m . We thereby conclude 
that, 

II, An optimum shape defined by the equations in (27) exists for a scaled length 
l > l*, when condition (28) holds and for l m >1 >1* when (28) does not. 

With a defined by (27), the coordinates of the points of tangency are defined by (17) and (18). 
The results are Tq( — 1 + Rq — Rx sin a, Ro cos a) and Tx(—Rx sin/3 + R\ sin a, h + Rx cos f 3 — 
Ri cos a). The scaled wave drag of the optimum forebody is 

*0,1,0) — — D 0 + Df + D n , (30) 

where D 0 -- ^![cosa(l — C °^ — ) — cos/?(l — C ° S — ], 

O o 

Df — sin 2 a[l + cos^ — ( Rq + Rx) cos a] 
and D n = Rq cos a[l — 

3 

Here Do,Df and D n denote the scaled wave drags of the optimum contour, the fillet and the 
nose. For hypersonic configurations, we have f3 << 1. With a < /3, we have D 0 = 0(/3 2 ), 
Df = 0(/3 2 ) and D n = |^ 0 + 0(Rq/3 2 ). If the nose radius is not very small, R 0 /3 2 , the 

wave drag of an optimum forebody is dominated by that of its nose, D„. and is only slightly 
less than that of the reference forebody by 0(/3 2 ). Since the scaled wave drag of a wedge 
of semi angle (3 ( Rq — 0) is sin 2 /3, the scaled wave drag of an optimum forebody will be 
greater than that of the wedge unless Rq is very small, 0(/3 2 ). For example, for a 4° half 
angle, /3 2 = 0.005 we have Ro/h ~ 0.005 which will be too small and most likely be ruled 
out due to aerodynamic heating. 

Consider two optimum forebodies having the same /3, R 0 , R x and h but different lengths 
l and l, as shown in Fig. 2. It is evident that the one which has the shorter length l has 
larger wave drag since 6(y) > 9(y). As the length l decreases the inclination of the common 
tangent a increases. When a reaches the upper bound /3, y x is equal to h and the optimum 
forebody coincides with the reference forebody of length l* with wave drag D*, which is 
independent of the fillet radius Rx- Again we note that the differences will be only 0((3 2 ) 
unless the scaled nose radius is very small. Thus we conclude: 


6 



III, The wave drag of an optimum forebody with R x > 0 is less than the wave 
drag D* of the reference forebody having the same (3, Rq and h while the length 
of the former is longer than the length l* of the latter. 

IV, If the nose radius is not very small, Rq ft 2 , the wave drag of an optimum 
forebody is dominated by that of its nose, D n and the reduction in the scaled wave 
drag for a longer forebody, l > l*, will be only 0(/3 2 ). 


4. AXI-SYMMETRIC OPTIMUM FOREBODIES 

In this section we seek the answers to the three items listed at the end of Sec. 3 for 
the axi- symmetric problems. In contrast to the two dimensional problems, we find that for 
each set of R, R x , Ro, (3 and l, with l lying within its upper and lower bounds, there are two 
optimum forebodies, called type o and b. The characteristics of these two types are studied 
in detail. For j = 1, Eq. (13) for the optimum contour becomes 

y sin 3 9 = hC for y 0 < y < yi, 
where hC = j/o sin 3 = yi sin 3 9 X . 

From (6), (8) and (32) we obtain an equation for 9 0 and 9 X . It is 

Ro cos 0 O sin 3 9 0 = C — (1 + R x cos/3 — R x cos 9 X ) sin 3 9 X . 

From (4), (5) and (6) we get the second equation, 

1+ Ro(l - sin0 o ) — -Ri(sin/3 — sin0i) 

= [1/A.] f cot 6 dy = [3/8]<7[flf(0i) -^(0o)], 

Jyo 

where di,®') = 2 esc 3 6 cot 9 — esc 9 cot 9 — log cot[0/2]. 

Since the optimum shape defined by (31) is convex while the fillet is concave, they can come 
in contact at only one point, if they do. That is Eqs. (33) and (34) can have at most one 
real root for 9 X with 

0 < 9 X < (3, while 0 < y x < h. (35) 

On the other hand both the spherical nose and the optimum contour (31) are convex. 
There can be more than one optimum solutions with different inclinations, 9o s. It is easier 
to see this by assigning Ro, R x ,f3 and 9 X and then solving for 9q from (33) and using (34) to 
define l. 


(31) 

(32) 

(33) 

(34) 
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From the equations in (43), we see that the scaled wave drag of the fillet, Dj, and the 
first term in the drag of the contour, Do are at most 0(/3 2 ). Using the optimum angle 9 X = 9 a 
from Eq. (48) for the optimum forebody of type a, we obtain from (43) the wave drag 

D a = 0(/5 2 ) + (R 2 / 4)[1 + 0(e 4 / 3 )] - (3/4 )R 2 0 e 2 ' 3 - 5 2 /4 ~ D n . (50) 

This says, that so long as the nose radius is not too small, say the leading contribution 

to the wave drag comes from the nose which is almost a hemisphere because 9 a e 2 / 3 « 1. 

For the optimum shape of type b we use Eq. (49) to define 9 X = 9b and obtain from (43) 
the wave drag 

D b = D } + (3/4 )y\ sin 2 9 X + 0(/? 6 ) = 0(/3 2 ). (51) 

This says that the contribution of the nose to the wave drag is of higher order 0(/3 6 ) relative 
to the total wave drag because the azimuthal angle of the spherical nose is 7r/2 — 9b ~ e. 

We used (43) and (44) to get (45), i.e., D a > Db. This inequality follows also from the 
explicit formulas (50) and (51). We note that we need the condition of Ro being not too 
small, (46), or rather Ro » ft, to establish (43) and (45). However, (43) - (45) hold for 
larger Rq. When, Rq = 0(1), the wave drag Db relative to D a is 0(/3 2 ). 

We conclude that • 

V, There are two optimum forebodies having the same l,h,R x ,&ndRo: type a, 
having an almost hemispherical nose, 9 a « 1 and type b, having a spherical 
nose with a small azimuthal angle, 7t/2 — 9b << 1. The wave drag of type b is 
the true minimum. 

VI, The wave drag of an optimum forebody of type b can be smaller the that of 
a conical body even when the nose radius R 0 = 0(1). The wave drag of type a 
will be greater than that of the conical body so long as the nose radius is not too 
small, i.e, Ro » (3. 


Since the two forebodies, a and b, have the same nose radius Ro, they have the same 
characteristics of aerodynamic heating near the stagnation point if the effect of entropy or 
vorticity gradient is negligible. For the forebody of type b, the vorticity gradient will be 
larger than that of forebody a and thus induce at larger heat flux [2,3]. The forebody b has a 
relatively smaller radius of curvature r b < Ro at the junction with its nose which could effect 
the boundary layer. Also the nose area could be too small such that the total aerodynamic 
fieating of the boundary layer near the stagnation point will not be sufficient to control the 
boundary layer downstream effectively. 

A realistic comparison of the optimum forebodies a and b should include both the wave 
drag and viscous drag. Additional theoretical and experimental investigations on the viscous 
drag, in particular on boundary layer control, are needed. 
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Now we shall find the lower and upper bounds for the scaled length J of the forebody for 
a given set of Ro, Ri, and/5. For the lower bound of l m (Ro, R\,0), we note that the maximum 
value of the constant C for the optimum contour is 

C m = min{sin A 0\ 3\/3i*o/16}. (52) 

When Cm = sin 4 /? we have 9 x — 0, and y x = h for both types a and b and for all i2i > 0. 
From C m we compute the two roots 9 a 9 b of Eq. (33) and the corresponding minimum length 
(L)a and (lm)b from Eq. (34). 

When C m = 3\/3^o/16 we have 9 a = 9 b = 7r/3 and 9 X defined by Eq. (33). Forebody a 
and b coincide with each other with l m given by Eq. (34). 

There is no upper bound on the length of an forebody of type a or b, if h + Ri cos 0 > 

Ro + R\ ■ 

In case that the above condition is not valid, there exists a least upper bound on the 
scaled length of an optimum forebody of type a. For this type, the circular fillet and the 
circular nose touch the optimum contour (31) on the opposite sides of the contour in a 
meridian plane. The scaled length reaches its maximum, {l u )a, when the circular nose is 
tangent to the circle of the fillet with y 0 = y x and the length of the optimum contour (31) 
vanishes. The condition on Ro, R x and 0 is 

\C 0 C 1 \/h = R 0 + R 1 . (53) 

with | C 0 Ci | defined by (27). The maximum (T\i)a is the same function of Ro,Ri and 0, 
given by (29) for the two dimensional problem. 

For the optimum forebody b, the circular fillet and the nose touch the optimum contour 
on the same side. There will be a least upper bound of ( l b when the circle of fillet intersects 
the x-axis. It occurs if h + R x cos 0 — R x <0 and the x-intercept is — i?x(sin 0 + sin 9 \ ) where 
9 X = cos -1 [(I + R x cos/3)]/i2i]. The least upper bound is 

(L) b = i?i(sin/3 - sin 9 X ). (54) 

This upper bound is not attainable by an optimum forebody b, since the optimum solution 
with J=Jb has a conical nose of cone angle 9 X with = 0. 

Using the geometry of the optimum contour (31) and that of the circular nose and fillet, 
we can show that the scaled wave drag of a forebody a or b, with given 0 and Ro increases 
as R x increase and/or l decreases. In the limit of R x — > oo, l approaches l* of the reference 
forebody. Thus we can conclude 

VII, The wave drag of an optimum forebody with a concave fillet, R\ > 0, for 
either type a or b, is less than the wave drag D* of the corresponding reference 
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forebody composed of a conical surface and spherical nose, with the same (3, Rq 

and h, while the length of the former is greater than the length l* of the latter. * 

5. CONCLUDING REMARKS _ 

Analytical solution for optimum shape of a two dimensional or axi-symmetric forebody « 

attached to a symmetric wedge or cone of half angle j3 is presented. The forebody is required = 

to have a circular nose of finite radius Rq defined by the allowable aerodynamic heating at the 1 

stagnation point. The forebody is joined to the truncated wedge or cone at a semi-thickness 
or base radius h through a circular fillet to avoid a discontinuity in slope. The radius of \ 

the fillet 72 x is specified with due consideration to avoid flow separation or shock boundary 1 

layer interaction.. Using the modified Newtonian theory, optimum forebodies are defined for 
given set of Ro,Ri, h,(3 and the length of the forebody l. f 

Using a reference forebody, composed of the extension of the truncated wedge or cone to I 

match with a circular or spherical nose of the same radius Rq, we show that the wave drag of 
an optimum forebody with a concave fillet, R\ > 0, ( a convex fillet R\ < 0) is less (greater) ■ 

than the wave drag D* of the reference forebody while the length of the former is greater = 

(less) that the length l* of the latter, see statement I. “ 

The optimum two dimensional and axi-symmetric forebodies with concave fillet junctions 
axe analyzed in detail in Secs. 3 and 4 respectively. The least upper and greatest lower 
bounds of the scaled length l of the forebody as functions of the parameters Ro,Ri,/3 with " 

h as the length scale are determined. The results are discussed under the realistic condition 
of /? <C 1. 

For the two dimensional cases, there is an unique optimum solution for a length 1 within 
the lower and upper bounds. The optimum forebody is defined by Eqs. (25) and (27) and 
its wave drag by Eq. (30). The wave drag is dominated by the part on the nose which 
is almost a semi-circle of radius R 0 so long as R 0 is not extremely small, i.e., R » (3 2 
and the reduction in drag relative to the drag of the reference forebody is small, 0((3 2 ) See 
statements II - IV. 

For the axi-symmetric cases, there are two optimum solutions of type a and type b for 
an scaled length l within the upper and lower bounds. The optimum forebody is defined by * 

Eqs. (31-34) and (40). The root of the last equation is given by (48) for type a and by (49) 
for type b. The wave drag is given by Eq. (43). The optimum forebody of type a is similar 
to the two dimensional case, in the sense that its nose is almost a hemisphere and its wave 
drag is dominated by that of the nose so long a,s Rq >> f3, The optimum forebody of type b 
has a small nose in the form of a spherical cap of small azimuthal angle while the radius Rq 
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can be 0(1). The scaled wave drag of the forebody of type b is 0(/3 2 ) and is dominated by 
that of the optimum contour and the fillet while that of the nose is of much higher order. 
The wave drag of type b can be smaller than the drag of the conical forebody of cone angle 
(3 even when Rq = 0(1). See statements V - VII. 

We note that our optimum solutions are based on the modified Newtonian theory. Since 
the surface pressure given by the modified Newtonian theory is in good agreement with the 
experimental data and numerical solutions of the full system of governing equations, we can 
introduce a small positive parameter 

8 « 1 

to denote the order of magnitude of the error of the surface pressure coefficient predicted 
by the modified Newtonian theory. Then the wave drag given by the experiment or by the 
numerical solution using our optimum forebody will differ from the true minimum by 0(6^). 
In other words, the “true” optimum shape will be small perturbation from our optimum 
shape. Instead of search for the “true” optimum shape directly by numerical solution it is 
much easier to determine the small perturbation numerically. 

It should also be noted that in addition to the wave drag we should also include the 
viscous drag. The latter requires a careful study of the boundary layer, flow transition 
and/or separation. Our optimum forebodies are simple and could serve as basic models for 
systematic experimental studies of the effects of the scaled length and nose radius on the 
total drag and their effects on the strong interaction of boundary layer with entropy layer 
and on the boundary layer control by utilizing aerodynamic heating near the stagnation 
point. 
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Figure 3. Optimum two dimensional forebodies with concave fillets. 
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